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^""^ Abstract 

C^ , In this paper we begin on a new lattice formulation of the non-linear change 

of variables called the Cho-Faddeev-Niemi decomposition in SU(2) Yang-Mills 
theory. This is a compact lattice formulation improving the non-compact lat- 

,^ , tice formulation proposed in our previous paper. Based on this formulation, 

we propose a new gauge-invariant definition of the magnetic monopole current 
which guarantees the magnetic charge quantization and reproduces the conven- 
tional magnetic-current density obtained in the Abelian projection based on the 

5^ I DeGrand-Toussaint method. Finally, we demonstrate the magnetic monopole 

dominance in the string tension in SU(2) Yang-Mills theory on a lattice. Our 
formulation enables one to reproduce in the gauge-invariant way remarkable 
results obtained so far only in the Maximally Abelian gauge. 
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1 Introduction 

In the previous paper [1], we have demonstrated that the gauge- invariant magnetic 
monopole can be constructed in the pure Yang-Mills theory without any fundamental 
scalar field. The success is achieved based on a new viewpoint proposed by three 
of us p] for the non- linear change of variables (NLCV), the so-called Cho-Faddeev- 
Niemi (CFN) decomposition which was first made by Cho [3] and has recently been 
readdressed by Faddeev and Niemi [1], see also [S]. We have for the first time formu- 
lated the NLCV on a lattice to perform non-perturbative investigations. Moreover, we 
have proposed a lattice construction of the gauge-invariant magnetic monopole which 
we called the lattice CFN monopole [1] . It was shown that the magnetic-current den- 
sity defined from the lattice CFN monopole reproduces the conventional monopole 
based on the DeGrand and Toussaint (DT) method [6] on a lattice. However, the 
previous construction of the lattice CFN decomposition and the resulting lattice CFN 
monopole have the following disadvantages. 

1. Gauge invariance on a lattice is manifest only to the first order of the lattice 
spacing e. 

2. The lattice artifact sets in at first order of the lattice spacing, since the finite 
difference is used in the calculations of differentiation. 

3. It is not obvious whether the resulting magnetic charge is quantized, since the 
monopole current is a real- valued variable. 

In this paper, we begin on a new lattice formulation of the NLCV or CFN decom- 
position, which resolves simultaneously all the issues raised above. This is nothing 
but a compact formulation of the Yang-Mills theory in terms of new lattice variables. 
Then we introduce a new definition of gauge-invariant magnetic monopole on a lat- 
tice based on the new formulation. The new formulation resolves simultaneously the 
first two issues, because we can define the electric and magnetic fields by a plaquette 
variable constructed from the new link variable as an element of a compact gauge 
group SU(2). The third issue is resolved and quantization of the magnetic charge is 
guaranteed, if the magnetic monopole is defined through the gauge-invariant fiux ob- 
tained by a plaquette variable. Moreover, we confirm dominance of the new monopole 
in the string tension, while it was first shown in [7] in the conventional Maximally 
Abelian (MA) gauge [8]. 

2 Lattice CFN variables or NLCV on a lattice 

In this paper we propose a natural and useful formulation on a lattice of the non-linear 
change of variables (NLCV) corresponding to the CFN decomposition [3],ll] in the 
continuum formulation. It is a minimum requirement that such a lattice formulation 
must reproduce the continuum counterparts in the naive continuum limit. In this 
stage, therefore, it is instructive to recall how the NLCV are defined in the continuum 
formulation. We restrict the following argument to SU(2) gauge group, since we wish 
to avoid technical difficulties and we give numerical results only for SU(2)1^ 



^ However, the following argument for SU(2) can be extended to the SU(3) case without many 
difficulties as will be shown in a subsequent paper. 
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In the continuum formulation [21[2], a color vector field n{x) = {nA_{x)) {A = 
1,2,3) is introduced as a three-dimensional unit vector field. In what follows, we 
use the boldface to express the Lie-algebra SM(2)-valued field, e.g., n(a;) := nA{x)TA, 
Ta = \(Ta with Pauli matrices u^ {A = 1,2,3). Then the sM(2)-valued gluon field 
(gauge potential) A^(x) is decomposed into two parts: 



A^(a;)=V^(a;)+X^(x), (1) 

in such a way that the color vector field n(x) is covariantly constant in the background 
field V^(x): 

= ^^[V]n(x) := d^x) - ig[V^ix), n{x)], (2) 

and that the remaining field X^(x) is perpendicular to n(x): 

n{x) ■ X^{x) = 2tr(n(a;)X^(x)) = 0. (3) 

Here we have adopted the normalization tr(T4Te) = ^Sab- Both n(x) and A^(x) are 
Hermitian fields. This is also the case for V^(a;) and X^(a;). 

By solving the defining equation ([2]), the V^(a;) field is obtained in the form: 

V,{x)=vUx)+V-tix) = c,ixMx)-tg-'[dM^),nix)], (4) 



where the second term V^(x) := —ig ^[d^n{x),n{x)] = g ^{d^n{x) x n(x))yiTyi is 
perpendicular to n(x), i.e., n{x) ■ V^{x) = 2tr(n(x)V^(x)) = 0. Here it should be 
remarked that the parallel part V|[(x) = c^(x)n(x), c^(x) = tr(n(x)A^(a;)) propor- 
tional to n(x) can not be determined uniquely from the defining equation ([2]). 

On a lattice, on the other hand, we introduce the site variable n^ constructed 
according to [1], in addition to the original link variable Ux,ij, which is related to the 
gauge potential A^(x) in a naive way: a 

f4,^ = exp{-iegA^{x)), (5) 

where n^ is Hermitian, nj_. = n^, and [/^^^ is unitary, Ul = U~]^. The link variable 
Ux,^ and the site variable n^ transform under the gauge transformation II [2] as 

Suppose we have obtained a link variable Vx^/j, as a group element of G = SU{2), 
which is related to the SM(2)-valued background field V^(x) through 



Vx,^, = exp{-iegYf,{x)), (7) 

must be unitary V^^^ = V~^. 



where V^(x) is to be identified with the continuum CFN variable (jll) and hence Vx^i^ 



^ In general, the argument of the exponential in ^ is the line integral of a gauge potential along 
a link from x to x + fi. We adopt this convention in agreement with our previous paper [l]. Note 
also that we define a color vector field n(a;) := nA{x)TA in the continuum, while n.^ :— u^cta on the 
lattice for convenience. 



A lattice version of Q and (jS]) is respectively given by 



'^x''x,^ ''x,^'^x+fii 



and 



tiin^U^^X,,) = 0- (9) 

Both conditions must be imposed to determine V^^^ for a given set of n^ and Ux^^- 
A lattice version of the defining equation ([2]) needs a lattice covariant derivative 
for an adjoint field. We adopt a definition of the covariant derivative for arbitrary 
background V^(x): 

D(^)[V]n, := e"i[K,;.n,+M - n.^J, (10) 

by the following reasons, i) When V"^^^ = 1, the derivative (TTO!) reduces to the 
(forward) lattice derivative ^jf^ria, := e"^[n^+^— n^;]. ii) The derivative (ITIII) reproduces 
correctly the continuum covariant derivative for the adjoint field up to 0{e): 

=e"^[n^.+^ - n^.] - ig[V ^{x)n.^+^ - n^V^(a;)] + 0(e), 
=dl:^n,-ig[V,{x),n,] + 0{e), (11) 

where we have used the O(e^) ambiguity in the last stepo iii) The derivative flTUj) 
obeys the correct transformation property, i.e., the adjoint rotation on a lattice: 

D(^)[V]n,. ^ n,{Dl^\V]nM+,, (12) 

provided that the link variable V^^^i transforms in the same way as the original link 
variable t/x./^: 

V.,, ^ n,v,,,nU^ = v;;^. (13) 

This is required from the transformation property of the continuum variable V^(x), 
see [2]. Therefore, we obtain the desired condition ([8]) between n^ and Vx^fj,. The 
defining equation (jH]) for the link variable V^^^ is form-invariant under the gauge 
transformation II, i.e., <1/'_^ = 14',^<+^. 

A lattice version of the orthogonality equation (^ is given by tr(nj.X^(x)) = or 

tr(n,exp{-ze(?X^(x)}) = tr(n,{l - ?e(?X^(x)}) + ^(e^) = + ^(e^). (14) 

This implies that the trace vanishes up to first order of e apart from the second order 
term. Remembering the relation X^(x) = A^(x) — V^(x), we can rewrite (TT^ into (Q 
in terms of n^ and U^^^- Note that the orthogonality condition (Q is gauge invariant. 
First, we proceed to solve the defining equation ([8]) for the link variable V^^^ 
and express it in terms of the site variable n^; and the original link variable Ux,n = 



^ Adopting another form instead of pip using the ambiguity of O(e^), — Ux+fi — rij. — 
iegV fj,{x)nx + inx+f_iegV^,{x) + 0{e^), which is rewritten as [1 + iegY ^{x) + 0{e^)\Ta.x = na;+^[l + 
ie,gV^(a;) + C(e^)], we obtain a relation V^^^n^ = n^+fiV^^^. However, this is nothing but the Hermi- 
tian conjugate of ([8]) and does not lead to a new condition. The definition of the covariant derivative 
could be improved e.g., by using a symmetric difference, as will be discussed in a separate paper. 



exp{—ieg Afj^{x)) , just as the continuum variable V^(x) is expressed in terms of n{x) 
and A^[x) in (jl]). The equation ([H]) is hnear in V^^^. Therefore, the normahzation of 
Vx^^ can not be determined by this equation alone. In general, unitarity is not guaran- 
teed for the general solution of the defining equation and hence a unitarity condition 
must be imposed afterwards. Moreover, the equation ([8]) is a matrix equation and it 
is rather difficult to obtain the general solution. Therefore, we adopt an ansatz (up 
to quadratic in n): 

Kr,M = U^,t^ + «n^f4,M + /?f^x,/.n^+/x + l^^xUx,^n^+^„ (15) 

which enjoys the correct transformation property, the adjoint rotation (TT^ . It turns 
out that this ansatz satisfy the defining equation^ if and only if the numerical coef- 
ficients a, (3 and 7 are chosen to be 

7 = 1, a = (3, (16) 

where we have used n^^n^ = n^n^aAO'B = '^^'^f ((^Asl + i^ABcO'c) = n ■ nl = 1. 



Second, substituting the ansatz ( 1151) with a still undetermined parameter a into 
the left-hand side of ([9]), we obtain 

tr(n,f/.,^V;t^) = a*[tr(l) + tr(n,.f/,,^n,.+^t/t^^)] = 4«* + 0{e'), (17) 

since 

tr(n,.t/,,^n,.+^f/i^) = tr{n.(l - zegAl^^){n., + e9j)n,.)(l + tegA^'^) + 0{e')} 
=tr{n,n, + ze^^n.n.A^f) - ze^^n^A^f^n, + en^dj^^n^ + 0(6^)} = tr(l) + ^(e^) (18) 

where we have used tT{nxdj^^nx) = + 0{e). Hence, the condition (^ leads to a = 
+ O(e^). Note that tr^n^+^V^^^Ux^fj) = gives the same condition as (^ by virtue of 
(E]). [If we imposed invariance on a lattice under the discrete global transformation, 
n^- -^ —iix respected in the continuum theory, we would have obtained a = = j3.] 
Thus we have determined \4,^ up to an overall normalization constant U 

Vx,f, = K,^[f/, n] = f4,^ + n^U^^^n^+f,. (19) 

Even in this stage, V^^^ is not necessarily unitary, since neither V^^^V^^^ nor V^^^V^^^ 
become a unit matrix 1. In fact, we have 

Vx,t,Vl^ = 21 + nMx,,,nx+iiUlf, + f/x-.^nx+Af^i^nx, 

VLVx,f. = 21 + t/J^^nxf/x,^nx+;. + n.+^Ul^n^U^^f,. (20) 

For Vx^fj, to become unitary merely by a normalization, the right-hand sides of Eqs. ( 120]) 
must be proportional to a unit matrix. We can show that this is indeed the case. |j 



"^ Uniqueness of the solution should be discussed separately. 

^ This special form has already been invented in a different context in the paper [9] in order to 
give the gauge-invariant lattice definition of Nambu magnetic monopolc with quantized magnetic 
charge in the SU(2) Higgs model on a lattice, although we have given more general scheme to find 
such a form in this paper. 

^ The parameterization, n^ — nl.ai and Ux^^j. = w".^ + Jw^^^fj, leads to nxUx,fj,nx+iiU^^^ = 
AqI + iAjGj with real variables Aq and Aj {j ~ 1,2,3): Aq = ('^x"x,/i)('^i+/i"i,;i) + 
{^l,i,nl - e'^''">ip)(w°,A."J+A' + e''™''?^i+^w"^)- Then, using a fact that C/^,^n^+/if/J;^^n^ 
and nxUx,^irix+fiU^ a ^'I'c Hcrmitian conjugate to each other, we find rixUx.^rix+fiU^ n + 
Ux,,,nx+f,Ul^nx = 2Aol. For a = 0, therefore, we have Vx,f,Vl^ = 2(1 + Aa)l = Vl^Vx,^,. 



Thus the unitary hnk variable \4,At[f^)i^] is obtained after the normahzation: 



K.M = V^AU, n] := V,jJ^tT[Vl^V,A- (21) 

Moreover, it is shown in the Appendix that the naive continuum hmit e ^ of the 
hnk variable fl?I]) reduces to the continuum expression (jl]). Therefore, we can define 
the gauge- invariant flux, Qp[U,n], (plaquette variable) by 

e.,^4U, n] := e-2arg(tr{(l + n,)V,,,%+f,,,V^+,^^VlJ/tT{l)). (22) 

This definition (!22|) is SU(2) gauge invariant due to cyclicity of the trace and the 
transformation property of the link variable and the site variable. This does not 
change thanks to the defining equation, even if we insert the factor (l+nj,) at different 
corners of a plaquette, e.g., Qx,f,u[U, n] = e-2arg(tr{V;,^(l + n^.+f,)V^+fi,uV^+u,,,Vl^}) . 
Finally, we show quantization of the magnetic charge as follows. We have con- 
structed a color vector field according to the adjoint orbit representation [T]: 

n. = e.asQl e, G 5f/(2), (23) 

which yields 

K-,A« = Ux,^i + n^.f/^,^na;+^ = Qx[Ux,i, + (^3Ux,^lCr3]Ql+^,, U^^^ := QlU^^^Q^+fi- (24) 
By representing an SU(2) element f/^;^^ in terms of Euler angles and Pauli matrices: 
U.,, := QlU,,,Qx+^ = e-3X./2e^<.2e./2g..3^./2^ q^ ^ [0, vr), y.,, xi e [-tt, tt), (25) 
the link variable Vx^fj, for a link £ = (x, fi) has the representation: 

Vx,ij. = 2 COS —Vx^fj,, 'Ke,/^ = ©z ( g ^_i(^^+^^)/2 1 ©x+AJ) (26) 



with the normalization factor given by y^trKr^^lxV = JAcos^^ = 2cosy. Thus 
the gauge- invariant flux ( 122|) is rewritten in terms of a compact variable $^ := {(pi + 
X£)/2e [-7r,7r) 

e,,^4f/,n] := e-^arg(tr{(l + n,)V;,MV;+A,.V;V.,;.V;tj/tr(l)) 

= ^-Vg(tr{(l + a3)(^^^;^^*^ ,-^eL^.)}^^(1)) 

= e-^argexpji'l'p} = [$p]mod2^, (27) 

where 

^P := {d<^)p = ^ $^ = <l>^_^ + '^x+t.,u - '^x+u,^^ - '^x,u- (28) 

It is important to remark that Qx,iJ,u on a lattice is a compact variable whose range is 
[— TT, tt), although it reduces to the continuum counterpart [I] which is non-compact 
variable taking the value (— cxd, oo) in the continuum limit. This fact is crucial to 
quantization of magnetic charge. In the unitary gauge, n^: = o"3, which corresponds 
to G^ = 1 in the above argument, Qx,tJ.u[U, n] agrees with 6x^f_iu in the DT field strength 
where the Abelian tensor 9^u{s) G [— 47r, 47r) C M is decomposed into the field strength 
part 0^u{s) G [— vr, vr) C M and Dirac string part n^^{s) G {—2,-1,0,1,2} C Z: 
(^nv{s) = 6^u{s) — 2Tm^y{s). It is known [13] that the elementary monopole defined in 
this way takes an integer-value [—2,2], since the Bianchi identity holds for 9^y{s). 



3 Numerical simulations 

As the algorithm for numerical simulations has already been given in a previous pa- 
per, we give a bit different explanation in this paper without repeating it. First of all, 
we generate the configurations of SU(2) link variables {Ux,^}, Ux,^j. = exp[— z5feA^(a;)], 
using the standard Wilson action based on the heat bath method. Next, we construct 
the color field variable n^. according to the following method. We minimize simul- 
taneously the two functionals F^mag and Fug written in terms of the gauge (link) 
variable U^,^ and the color (site) variable n^,: 

F„MAG[f/,n;fi,e] :=5]tr(l- Vf/x.^V+M^'t^U' (29) 

Fllg[U-M :=Etr(l-^t/x,;.), (30) 

with respect to two gauge transformations: ^Ux,fi '■= ^xUx,^^x+ii for the link vari- 
able Ux,^ and ®n^. := 6a;n^°''Gj, for an initial site variable n^^-* (we can choose the 
initial value n^'') = a^) where the gauge group elements Q^ and 6^; are indepen- 
dent SU(2) matrices on a site x. Then we can determine the configurations ® n^ 
and ^"Ux,^ realizing the minimum of the first functional: min^ e F^mag]!^-, ^i ^' ®] ~ 
F„A'/yiG[f^) n; r2*, 6*], up to a common SU(2) transformation G^- This is because 
the "common" gauge transformation G^ for 6* and ^2* does not change the value of 
the functional FnUAGp, n; Vt, 6], i.e., FnUAGp, n; VL\ 9*] = FnMAGp, n; GVL\ GG*], 
since tr(~n2,. Ux,^~'iix+^i Ul^^) = tT{Gx~'nxGx+fj, ■ Gx+ii Ux,ij.Gx+^ ■ Gx+ij.~'n.x+tiGx+^ ■ 
Gx+f/^U^^^Gl). This degrees of freedom for the SU(2) gauge transformation are fixed 
by minimizing the second functional Fllg[U', ^] such that the configuration ^ f/^_^ re- 
alizes the minimum of the second functional: min^ Fllg[U; fi] = Fllg[U] fi**]. Thus, 
imposing simultaneously two minimizing conditions removes the SU(2) ambiguity and 
the color field configuration are decided as n^; := Bajii^P^Oj. with Qx = ^Ti^D^^^l 

with Gx = ftT{n*x)-\ 

Once the configurations of the color vector field {n^;} are generated according to 
the method explained above together with the configurations of SU(2) link variables 
{Ux.ii}, we can construct {\4,/^[f^, n]} from (12T1) . The numerical simulations are per- 
formed on an 8^ lattice at p =2.2, 2.3, 2.35, 2.4, 2.45, 2.5, 2.6 and on 16^ lattice at 
/5 = 2.4 by thermalizing 3000 sweeps respectively. 

3.1 Quantization of magnetic charges 

We construct the gauge-invariant field strength (1^ to extract configurations of the 
magnetic monopole current {kx^^} defined by 

47r 47r 

This definition agrees with our definition of the magnetic monopole in the continuum 
(divided by 2tt). This definition of the monopole current should be compared with 
the conventional monopole current on a lattice defined according to DeGrand and 
Toussaint [6] through link variables on the dual lattice [1]: 



kli[-^) — n^p.vpcjt^v'npcfyS -\- H) — ^pupa(^u^pcj\,S ~\~ l^j- 
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The monopole current k^{s) defined in this way becomes an integer- valued variable, 
since integer- valued variables np„ are used to count the number of Dirac strings going 
out through a plaquette. 

In our formulation, on the other hand, we have only the real variable Bp[?7, n] 
at hand, and we are to calculate the monopole current using the final term in (13T1) . 
Therefore, it is not so trivial to obtain the integer- valued kf^{s) from the real- valued 
0p[?7, n]. To check quantization of the magnetic charge, we have made a histogram 
oi K{s,fi) := 2nkfj_{s) = ^e^upadiyQpa{x + fi), i.e., magnetic charge distribution. Note 
that K{s, fi) should become a multiple of 27r if the magnetic charge is quantized. Our 
numerical results show that K{s,fi) is completely separated into or ±2tt within 
an error of 10~^°, see Table [H We have checked that the data in Table [T] exhaust 
in total all the configurations A^ = 4 x 8^ = 16384, because the number A'^; of links 
in the d-dimensional lattice with a side length L is given by Ni = dL'^. This result 
clearly shows that the magnetic charge defined anew is quantized as expected from 
the general argument. We have observed that the conservation law of the monopole 
current holds, since the number of +2tt configurations is the same as that of — 27r 
configurations. In contrast. Table [1] shows that quantization does not occur for the 
(old) CFN monopole constructed in the previous paper [T] using Hf^^, alone. 



Table 1: Histogram of the magnetic charge (value of K{s,fi)) distribution for new 
and old monopoles on 8^ lattice at /3 = 2.35. 



Charge Number (new) Number (old) 



-7.5—6.5 








-6.5~-5.5 


299 





-5.5~-4.5 





1 


-4.5~-3.5 





19 


-3.5—2.5 





52 


-2.5—1.5 





149 


-1.5—0.5 





1086 


-0.5-0.5 


15786 


13801 


0.5~1.5 





1035 


1.5~2.5 





173 


2.5-3.5 





52 


3.5~4.5 





16 


4.5~5.5 








5.5-6.5 


299 





6.5-7.5 









3.2 Magnetic-current density 

Storing 30 configurations of the monopole current, we have calculated the magnetic- 
current density Pmon defined by 



Pti 



-<j:\k,{x)\>/{W). 



(32) 



X,IM 



The results are summarized in the second column of Table [21 This should be compared 
with the first column of the conventional DT magnetic-current density [TT] and the 



third column of the old (naive) CFN one calculated from H^y alone [T]. The DT 
magnetic-current density is calculated from the photon potential A^ extracted by 
using the Abelian projection [12] in the MA gauge. The new monopole agrees with 
DT monopole in density to high accuracy better than the naive CFN monopole. 
The similarity of the magnetic charge and current-density distributions between the 
DT monopole and the new monopoles suggests the dominance of the new magnetic 
monopole in the string tension, as will be discussed shortly. 

Note that the new monopole always keeps gauge invariance by including both the 
electric and magnetic contributions coming from c^ and H^^ respectively. In contrast 
to the naive CFN monopole, however, the new monopole is difficult to be separated 
into the electric and magnetic parts in this sense. 

Table 2: j3 dependence of magnetic-current densities for DT monopole, new and old 
CFN monopole. 



/? 


DT monopole 


new 


old 


2.2 


0.0825(4) 


0.0821(7) 


0.0733(5) 


2.3 


0.0515(3) 


0.0514(12) 


0.0540(6) 


2.35 


0.0379(3) 


0.0371(7) 


0.0423(8) 


2.4 


0.0263(4) 


0.0257(7) 


0.0326(8) 


2.45 


0.0167(6) 


0.0172(9) 


0.0233(10) 


2.5 


0.0096(5) 


0.0108(8) 




2.6 


0.0040(2) 


0.0036(3) 





3.3 Monopole dominance of the string tension 

In order to study magnetic monopole dominance in the string tension, we proceed to 
estimate the magnetic monopole contribution (WmiC)) to the Wilson loop average 
{Wf{C)), i.e., the expectation value of the Wilson loop operator. We define the 
magnetic part Wm{C) of the Wilson loop operator Wf{C) as the contribution from 
the monopole current /^^(s) to the Wilson loop operator: 13 

Wm{C) =exp|27rz^A;^(s)iV^(s)|, (33) 

N,{s) =Y.^l\s - s')\e,^p,d^Sf,^{s' + fi), %Sl^{s) = J,{s), (34) 

s' 

where N^{s) is defined through the external electric source J^{s) which is used to 
calculate the static potential: d' denotes the backward lattice derivative d f{x) = 
f{x) — f{x — /i), S^^{s) denotes a surface bounded by the closed loop C on which 



"^ The Wilson loop operator Wf (C) is decomposed into the magnetic part Wm (C) and the electric 
part We{C), which is derived from the non- Abelian Stokes theorem, see Appendix B of [T3]. In this 
paper, we do not calculate the electric contribution {We(C)) where We{C) is expressed by the electric 
current j^ = d^Ff^^. 



1.4 












flJII 
monopole 


■ 


1.2 












T 


- 


1 
> 






I'-"" 


,^'-"f"" 


..^^'*""" 


'-'"-' 


- 


0.6 




.w' 










- 


0.4 


/' 








.— '^""" 




- 


0.2 




.-^'•'-" 


-^^-'--""" 








^ 


1 


2 


3 


4 


5 


6 7 


8 



Figure 1: The full SU(2) potential Vf[K) and the magnetic-monopole potential 
Kn(-R) as functions of i? at /5 = 2.4 on 16^ lattice. 

the electric source J^i^s) has its support, and A^^(s — s') is the Lattice Coulomb 
propagator. We obtain the string tension by evaluating the average of ( !33l) from the 
generated configurations of the monopoles {/i;^(s)}. 

The numerical simulations are performed on 8^ lattice at /? = 2.3, 2.4, and IG'^ 
lattice at /3 = 2.4 by using 50 configurations in each case with 100 iterations. In 
particular, we have used 100 configurations for the calculation of the full string tension 
at /3 = 2.4. The Wilson loop average (Vr(i?, T)) for a rectangular loop C = (-R, T) 
with side lengths R and T is calculated by varying R from 1 to 4 for a fixed T (T = 4) 
on 8^ lattice, and from 1 to 7 for a fixed T (T = 7) on 16"^ lattice. Then we have 
calculated the respective potential Vi{R) from the respective average {Wi{C)): 

V,{R) = -\og{{W,{R,T))/{Wi{R,T-l))} {i = f\m). (35) 

Fig. [1] shows the full SU(2) potential V/(i?) and the magnetic-monopole potential 
Vm{R) as functions oi RaX (3 = 2.4 on 16^ lattice. El Moreover, the numerical potential 
is x^ fitted to the form with a linear term, the Coulomb term and a constant term: 

Vi{R) = aiR-ajR + Ci, (36) 

where a is the string tension, a is the Coulomb coefficient, and c is the coefficient 
of the perimeter decay: {Wi{R, T)) ~ exp[— o"j-RT — Cj(i?, + T) + aiT / R +•■■]. The 
results are shown in Table El 

Table 3: String tension and Coulomb coefficient I 

13 CFf O/ (Jm Orn 

2.3(8^) 0.158(14) 0.226(44) 0.135(13) 0.009(36) 
2.4(8^) 0.065(13) 0.267(33) 0.040(12) 0.030(34) 
2.4(16^) 0.075(9) 0.23(2) 0.068(2) 0.001(5) 



Fig. [Hand Table [S] show that the magnetic-monopole potential Vm{R) has a dom- 
inant linear term and a negligibly small Coulomb term with small errors. Conse- 
quently, the string tension a^. is obtained within small errors. In order to obtain the 



^ It should be remarked that in Fig. 2 of Stack et al. [7], a constant term was added to a 
magnetic-monopole potential to compare the i?-dependence of the potential. 
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full SU(2) potential Vf{R), on the other hand, various techniques have been used to 
eliminate as much as possible the unphysical short-distance fluctuations [IH]. Espe- 
cially, we have used (APE's) smearing method [13] as a noise reduction technique. In 
this method, the Wilson loop average is calculated by using the block-link variable 
U'^ constructed from the link variable Ux,fj, according to the procedure: 

K, = U.,, + « E f/.,.^x+p,Mf/i+A,- (37) 

where a is a parameter to be adjusted later and the blocking is performed only for 
spacial directions, not for the time direction. Of course, the obtained U'^^^ is not 

unitary and hence it should be projected to be a unitary matrix t/^^^ by a standard 
method. This procedure is repeated TV times iteratively starting from initial conflg- 
urations obtained by the heat bath method. The parameters A^ and a are chosen so 
that the Wilson loop average takes the maximal value. In our calculations for the full 
SU(2) potential on 16^ lattice, we have used 100 conflgurations (twice the number of 
the other cases) to decrease statistical errors of af and aj, and adopted A^ = 15 and 
a = 0.2 according to [TTj. The gauge invariance is preserved by this procedure, since 
the block-link variable U'^ {U^ ) has the same transformation property as the link 
variable U^^fj,- 

In fact, we find that the monopole part am reproduces 85% of the full string 
tension o"/ at /3 = 2.3 on 8^ lattice and that am reproduces 91% of o"/ at /3 = 2.4 on 
16^ lattice. Thus, we have confirmed the magnetic monopole dominance in the string 
tension using our magnetic monopole in the gauge-invariant way. 

Table 4: String tension and Coulomb coefficient II (reproduced from [7]) 



(3 CTf O/ CrpTm aPTm 

2.4(16'^) 0.072(3) 0.28(2) 0.068(2) 0.01(1) 

2.45(16^) 0.049(1) 0.29(1) 0.051(1) 0.02(1) 

2.5(16^) 0.033(2) 0.29(1) 0.034(1) 0.01(1) 



For comparison, we have shown in Table H] the data of [7] which has discovered the 
monopole dominance for the first time on 16^ lattice where aoTm reproduces 95% of 
af. lI Here aoTm and aoTm denotes the conventional monopole contribution extracted 
from the diagonal potential A'i using Abelian projection in MAG. In particular, the 
comparison of the data on 16^ lattice at /5 = 2.4 between Table [Hand Table O reveals 
that the monopole contributions have exactly the same value between the conventional 
DT monopole and our magnetic monopole. This is because the monopole part does 
not include the Coulomb term and hence the potential is obtained to an accuracy 
better than the full potential, as pointed above. 



^ Stack ct al. [7j have have performed numerical simulations using 500 configurations with 20 
iterations on 16"* lattice in the calculation of the potential. This will lead to less statistical errors for 
the full string tension by a factor l/VS than ours, since we have used 100 configurations at /? = 2.4. 
Indeed, our data are consistent with this estimation. In addition to the smearing method, they 
used also the multihit method [H] to decrease noises in the calculation of the full string tension. 
Incidentally, they have measured the potential up to i? x T := 7 x 10 on 16^ lattice. 
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4 Conclusion and discussion 

In this paper, we have proposed a new formulation of the NLCV of Yang- Mills theory 
which was once called the CFN decomposition. This resolves all drawbacks of the 
previous formulation [1] on a lattice. This compact formulation enables us to guaran- 
tee the magnetic charge quantization. The new monopole dominance has been shown 
anew in the string tension. In this paper, the magnetic charge quantization and the 
magnetic monopole dominance in the string tension are confirmed in the gauge in- 
variant way, whereas they have been so far shown only in a special gauge fixing called 
MA gauge which breaks the color symmetry explicitly. 

A suitable definition of the link variable exp{—ieg^fj_{x)) for X^(x) will be given 
with relevant numerical results in a separate paper. Extending the promising formu- 
lation proposed here for SU(2) to SU(3) will be given in a subsequent paper. 
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A Naive continuum limit 

We consider the naive continuum limit e ^ of the quantities on a lattice defined 
above to see that they are good definitions on a lattice. 

A.l Continuum limit of 14^ 

The definitions of the gauge field [4^^ = e"*''^^''^^) (A^(a;) = A^(x)^) as a link 
variable and the adjoint scalar field n^ = n'^a"' as a site variable yield 

= nl + en^d^n^-iegn^Af,{x)n^ + 0{e^), (38) 



where we have used the relations: (o"")^ = 1 (no summation), {a", a} = 26 1 



[cr°, a"] = 2ie'""'a'', which follow from a^a" = S^"! + ie^'^'a". The first term of (EH]) 
reads n^ = n'^n'^a"'a^ = n"?T.^(5'^''l -|- ie°'^'^a'^) = 1, where we have used (n")^ = 
n°'n°' = 1, e'^'^^n'^n'' = 0. The second term of (138!) reads n^df^tix = n%dfj_n\a°'a^ = 
—i^o-bc ^g^^a-^^b^^c _. _2i(y(B^(x), where we have used n%d^j_n% = which follows from 
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differentiating (ra")^ = 1. Ttie tfiird term of flHHl) reads 

= -iegAl{x)ay2 + z6^KA;,(x))<a^ 

= -iegA^{x) + 2iegc^{x)ii{x), (39) 

wliere we liave used . twice and g'^^'^g'^ce _ ^ac^he _ ^ae^bc^ ^^^ defined c^{x) by c^{x) = 
?7,"A°(x). Substituting tliese results into (1HH1) . we obtain 

n^t/^,^n^+^ = 1 + ie^(A^(a;) - 2c^(x)n(x) - 2B^(a;)) + ^(e^) 

= l + tegiA^ix)-2V^ix)) + 0{e^). (40) 

Tfierefore, we arrive at 

= (1 - ^e(7A^(x)) + (1 + ieg{A^{x) - 2V^(x)) + ^(e^) 

= 2(1 - tegV^ix)) + 0{e') = 2K,^(f/, n). (41) 

Hence, V'^(x) in the continuum limit is calculated from egYfj,{x) = — ^(t4-,/i, — Vx,^l)■ 

A. 2 SU(2) gauge-invariant field strength Qx^fiu 

The plaquette variable Vp constructed from the link variable V^.,^ reads 

Vp := V.,>.Vx+i.,uVl+.,^Vl. = exp{ze^g^,4Y]}, (42) 

with the field strength ^^^[V] for V^ defined by ^^.^[V] = 9/,V^-9^V^-z^[V^, V^]. 
Then we have 

tr{n^V>} = tr(n,.) + «eVr(nx^M-[V]) + 0(e^) 

= ze^gtT{nx{d^y, - d,Y , - ig[Y,,\,])} + O(e^) 

= ie^gii{l){n ■ [d^{ci,n + g~^d„n x n) - d^{c^n + g'^d^n x n) 

+g{Cfji + g~^d^n X n) X {c^n + g^^duU x n)]} + O(e^) 
= iehi{l){d^c,-d,c^ + g-'n-{d^nxd,n)} + 0{e^), (43) 



where we have used n^. = 1, {a x h) x {a x c) = {a ■ h x c)a . On the other hand, (142|) 
leads to tr{Vp} = tr(l) + C(e^). Finally, we find that Qx,^ly has the same form as the 
't Hooft tensor [T6l: 



tr{(l + n^.)^p}/tr(l) = I + ie^{d^Cy - d^c^ + g'^n ■ d^^n x d^n)} 

= exp{ie^{df,Cu - d^c^ + g'^^n ■ d^n x dyu)} 
= explie^e^,^^. (44) 

Thus we have constructed a lattice version of the gauge-invariant field strength: 

©x,/.!. ^ df,c^ - d^c^ + g~^n ■ {d^n x d^n) = G'^i.(x). (45) 
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